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Stress  analysis  procedures  are  presented  for  fibre-reinforced  metallic 
rocket  motor  cases  using  finite-element  and  photoelastic  techniques. 

These  techniques  enable  analyses  of  the  initial  winding  process,  inertia  i 
and  preasurlsatlon  due  to  firing  to  be  performed.  Certain  simplifications 
and  assumptions  are  Introduced  to  resolve  the  complicating  effects  of  the 
fibre  overwind  and  the  in-flight  loadings.  A  comparison  between  the  two 
methods  for  a  typical  motor  is  considered. 
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1  INTRODUCTION 

Hie  selection  of  materials  with  high  specific  strengths,  le  the  ratio  of 
ultlaate  tensile  strength  to  specific  gravity,  Is  an  important  feature 
in  the  design  of  rocket  motor  cases.  Maraglng  steel  for  example,  with  an 
ultimate  tensile  strength  of  2200MPa  and  a  specific  gravity  of  8,  Is  one 
of  the  highest  specific  strength  materials  suitable  for  'single'  material 
designs.  Attempts  to  use  the  even  higher  strength  steels  have,  however, 
generated  serious  quality  assurance  problems  associated,  for  instance, 
with  their  reduced  ductility  and  some  motors  have  failed  In  a  brittle 
manner.. 

Alternatively,  with  the  choice  of  suitable  constituent  materials,  composite 
rocket  motor  cases  can  be  manufactured  with  specific  strengths  exceeding 
that  of  any  homogeneous  design.  A  lightweight,  load-bearing,  almlnlun 
tube,  circumferentially  reinforced  with  a  high-strength  polyamide  fibre, 
is  a  particular  example.  Furthermore,  if  the  fibre  overwind  Is  applied 
under  tension,  the  case  is  put  Into  compression.  This  defers  the  onset 
of  tensile  stresses  and  thereby  effectively  Increases  the  firing  pressure 
of  the  motor.  In  spite  of  these  advantages  prestralnlAg  has  potential 
drawbacks.  For  example,  experimental  tests  have  shown  that  the  strength 
and  modulus  of  the  fibre  are  reduced  with  prolonged  periods  of  loading. 
Nevertheless,  provided  that  adequate  allowances  are  made  for  such 
degradations  at  the  design  Stage,  the  considerable  Increase  In  the 
effective  strength  usually  outweighs  these  drawbacks. 

In  the  past  the  stress  analysis  of  flbre-relnforced  metallic  rocket  motor 
cases  has  been  performed  empirically  (Ref  1),  but  during  recent  years 
the  authors  have  been  Involved  in  the  development  of  theoretical  stressing 
procedures  (Refs  2-6).  This  work  has  involved  analyses  of  the  Initial 
winding  process  for  both  constant  and  variable  winding  strains,  temperature 
excursions  about  the  ambient  winding  temperature  and  the  response  of  the 
motor  when  subjected  to  Its  firing  pressure. 

For  the  thln-walled  cylludrlcel  portion  of  the  motor  analytlcel  stress 
solutions  have  been  developed,  whereas  for  the  thicker  and  geometrically 
more  complex  regions  of  the  motor  finite-element  procedures  have  been 
applied.  Fhotoelaatlc  stress  analysis  procedures  have  also  been  developed 
for  fibre-reinforced  rocket  motor  cases  (Refs  7-8)  so  that  the  three-*.. 
dimensional  effects  of  the  geometry,  which  cannot  be  readily  analysed,  *».*, 
using  the  theoretical  procedures,  can  be  studied.  In  this  report  the 
theoretical  and  photoelastlc  techniques  which  have  been  used  to  analyse 
fibre-reinforced  metallic  rocket  motor  cases  ere  reviewed  and  their 
respective  stress  solutions  compared  for  a  typical  motor. 

2  TH1N-HALLED  SECTION  ANALYSIS 

2.1  Overwinding  Analysis 

When  n  layers  of  fibre  overwind  ere  wound  clreiaeferentlally  at  an  Initial 
winding  strain  cw  on  a  metal  case  of  radius  R  the  case  contracts  radially 
owing  to  the  external  preaaura  of  the  overwind  and  la  thereby  put  In 
compression.  In  the  absence  of  an  Internal  pressure  the  axial  stress  Is 
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zero  and,  aaaualng  that  the  radial  contraction  of  the  motor  due  to  the 
application  of  each  fibre  layer  la  uniform  through  the  overwound  section. 
It  hae  been  shown  (Ref  5)  that  the  elastic  pre-stress  In  the  l'th  fibre 

layer  and  the  case,  and  *«*P«etlvely,  are  given  by 


oef(1'n)  -  Eft^tl-AFjU.n)]  ,  i-  1,2 . . 


(1) 


«0c(n>  *  “*cewAFl(A»n) 


where 


Ft(A,n) 


n 

z  „ L. 


(2) 


(3) 


with  the  stiffness  parameter  A  denoting 


In  these  equations  the  quantities  Ef  and  Ec  are  the  Young's  modulus  values 
for  the  fibres  and  case  respectively,  tf  Is  the  effective  thickness  of  o 
single  fibre  layer  and  tc  Is  the  wall  thickness  of  the  metallic  case. 

for  ease  of  computation  It  Is  convenient  to  approximate  the  simulation 
tern  (le  Fi(A,n))  appearing  In  Equations  1  and  2  by  an  analytical 
expression.  Using  the  Euler-Maclaurin  formula  (Ref  9)  and  neglecting  all 
first-order  terms  and  above.  Equation  3  modifies  to: 


Fj(A,n) 


I  ln[ii*S.]  +  1  +  1  .  . 

A  1 14A1 J  2(14An)  2(14Al!) 


(5) 


Detailed  computations  have  shown  that  for  typical  A  values  (le  A  1)  the 
errors  associated  with  this  approximation  are  very  small  (of  the  order 
0.05Z). 

Once  the  stresses  In  the  fibres  and  the  case  have  been  derived,  the 
corresponding  elastic  strains  follow  from  the  uniaxial  Nookean  stress- 
strain  relationship. 

2.2  Preseurlaatlon  Analysis 

When  the  rocket  motor  case  Is  subjected  to  Its  firing  pressure  It  Is  a 
requirement  that  the  motor  must  withstand  all  ln-fllght  loads  without 
failure  or  excessive  deformation.  Although  these  ln-fllght  loads  may 
contain  components  due  to  lnertls  and  aero-heat,  the  dominant  load  will 
always  ha  the  firing  praasure. 
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In  the  overwinding  analysis  the  axial  stress  was  zero.  However,  when  the 
aotor  is  subjected  to  Its  firing  pressure  the  axial  stress  will.  In 
general,  be  non-zero  and  la  given  by 


(R2-Re2> 

°*CP  "  P  2Rt  c  ’ 


(6) 


where  Ra  la  the  effective  Internal  radius  (Ref  3)  of  the  throat  Insert, 
taking  Into  account  the  non-uniform  pressure  distribution  acting  over  the 
nozzle  region.  When  an  allowance  for  Inertia  effects  la  Included,  Re  will 
contain  an  additional  component  due  to  the  axially  varying  Inertia  force. 

The  circumferential  atreasea  In  the  fibres  and  the  case  are  also  modified 
when  the  motor  la  subjected  to  its  firing  pressure.  It  can  be  shown  from 
a  previous  analysis  (Ref  3)  that  the  modified  elastic  clrcwferential 

stresses  In  the  i'th  fibre  layer  and  the  case,  cef  (1,n)  and  oecp^ 
respectively,  are  given  by 


°efp<1,n)  *  EfZjl-AF^A.n)] 


+  .  P"A . 

tf(l+nA) 


l,2p  e.etn# 


(7) 


«0cp(n)  -  -VwVA-"> 


.£*  .  [i+: 

1  imi  S  L 


tc(l-HiA) 


nvcA 


(8) 


where  vc  denotes  the  Poisson's  ratio  of  the  ZMtalllc  case  and  p  Is  the 
firing  presaura. 

As  previously,  once  the  stresaea  have  bean  obtained,  the  elastic  s '.rains 
follow  from  the  Hookean  atresa-atraln  relatione. 

3  TNI  CTC -WALLED  SECTION  ANALYSIS 

In  the  preceding  section  simple  analytical  stress  solutions  were  presented 
for  the  thln-walled  sections  of  the  meter.  For  the  thick -walled  end 
geometrically  complex  regions  of  the  ease  a  more  eophietlcated  stress 
analysis  approach  is  required. 

Por  complex  geometries  the  required  strees  analysis  can  be  performed 
either  theoretically  or  experimentally.  Tor  a  theoretical  analysis  flnlte- 
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element  technique*  have  to  he  adopted,  whereas  for  an  experimental  study, 
strain  gauging  or  photoelastic  procedures  can  be  applied.  For  experimental 
studies  the  photoelastic  technique  is  usually  preferred  since  It  provides 
stress  Information  at  both  the  surface  and  Interior  points  of  the  structure. 
It  also  has  the  advantage  of  providing  detailed  Information  of  the  stress 
variations  In  the  regions  of  high  stress  concentrations  taking  Into 
account  full  three-dimensional  effects  of  the  geometry. 

The  finite-element  and  photoelastic  techniques  are  well  documented  (kefs 
10  and  11)  and  are  both  capable  of  providing  accurate  stress  predictions. 
However,  the  aenner  in  which  these  techniques  ere  applied  to  the  overwound 
rocket  motor  case  Is  not  entirely  straightforward  due  to  the  complicating 
effects  of  the  fibre  overwind  and  the  ln-f light  loadings.  Both  methods 
of  analysis  are  now  discussed. 

3.1  Finite-Element  Analysis 


For  a  rigorous  stress  analysis  of  the  flbre-relnforced  metallic  rocket 
motor  case  it  Is  necessary  to  employ  specialised  finite-element  programs 
which  provide  facilities  for  accurately  modelling  the  reinforcing  properties 
of  the  fibre  overwind  and  the  pre-stress  produced  during  winding.  Although 
such  programs  are  under  -'evelopaient,  they  are  not,  as  yet,  available  for 
use  and  It  has  been  necessary  to  develop  an  alternative,  albeit  approximate, 
aolutlon  procedure  using  existing  finite-element  programs. 

In  the  theory  of  finite-element  stress  analysis  there  exists  a  technique 
known  as  sub-structuring  which  enables  the  analysis  of  complex  structures 
to  be  simplified.  This  technique  Involves  the  replscement  of  part  of  a 
structure  by  a  set  of  equivalent  forces,  or  displscements,  applied  to  the 
Interface  boundary  of  the  retained  portion  (Ref  10).  For  the  overwound 
motor  It  follows  therefore  that  the  overwind  can  be  removed  and  structurally 
represented  by  an  external  pressure  equal  In  magnitude  to  the  Interface 
pressure  between  the  first  fibre  layer  and  tha  metal  case  acting  over  the 
region  previously  occupied  by  the  fibre  overwind.  The  remaining  structure, 
a  homogeneous  Isotropic  matal  case,  can  then  be  analysed  by  any  conventional 
axlsymmetrlc  finite-element  progrsn. 

At  this  stag*  it  would  seem  logical  to  evaluate  the  flbre/cas*  Interface 
pressure,  py,  by  first  calculating  the  circumferential  streas,  qj,.'11', 
from  Equation  8  and  then  using  the  thln-wsllod  cylinder  equation  to  give 


h 


P 


*c 


(8) 


However,  la  sltmatloas  where  the  Interface  preeeure  py  and  the  applied 
internal  preeewre  p  are  alaeet  equal.  It  has  been  feeed  from  am  error 
analysis  for  the  drew* sraatlal  stress  (see  dames  A)  that  as  py  approaches 
p  a  aaell  error  la  the  determination  of  py  will  load  to  a  large  error  In 
the  circumferential  s truce  calculation.  Although  the  above  equation 
yields  aa  Interface  pressure  within  a  few  par  cent  of  that  predicted  by 
the  enact  Lama  equations  (Ref  Id),  serious  errors  wars  observed  la  the 
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finite-element  stress  solutions  (in  excess  of  twenty-five  per  cent)  when 
this  estimate  of  p^  was  used.  To  ensure  that  accurate  results  are 
generated  it  is  vital  that  the  more  precise  relationship  based  on  the 
Lame  equations  (see  Annex  A) 


**b 


r  2Ri2 

vc(R12-Re2 

LEc(Rbi-R1'i) 

Ec(Rb2-Riz  > 

r  Rv,2+r<2 

+  RrMh2  ] 

1  VV-*i‘£) 

Ef(V-Rb'J)J 

A  EcewAF1(A,n)tc 
_ 


(10) 


is  used  in  the  finite-element  analysis,  where  Rj  and  Rb  denote  respectively 
the  internal  and  interfere  radii  of  the  metallic  case  and  R0  is  the  external 
radius  of  the  overwound  motor,  ie  Rb  +  ntf. 

In  this  approximate  analysis  the  overwound  section  of  the  motor  is  divided 
into  a  series  of  cylinders,  the  lengths  of  which  are  selected  to  take 
into  account  the  non-uniform  axial  stress  variation  due  to  inertia  and 
possible  variations  in  the  total  thickness  of  the  fibre  overwind  and  the 
case.  For  each  cylinder  the  quantity  Re  appearing  in  Equation  10  is 
calculated  from  Equation  6  in  terms  of  the  axial  stress.  This  axial 
stress  consists  of  two  components,  one  which  is  proportional  to  the 
constant  axially  resolved  Integrated  force  due  to  the  non-uniform  internal 
pressure  acting  over  the  nozzle  region,  and  one  due  to  inertia,  which  is 
axially  non-uniform.  It  should  be  noted  that  in  the  vicinity  of  the  ends 
of  the  overwound  section  of  the  case  the  cylinder  approximation  is  not 
strictly  valid,  but  experimental  results  Indicate  that  this  method  of 
analysis  provides  a  reasonable  working  model  for  this  theoretical  study 
(see  Section  4). 

Since  Rg  is  derived  from  the  axial  stress  distribution  it  follows  that 
the  fibre/case  Interface  pressure  (see  Eqn  10)  is  a  function  of  the  stress 
state  being  calculated.  The  solution  procsss  is  therefore  iterative, 
requiring  an  initial  approximation  for  Re,  and  hence  the  Interface 
presaure,  for  each  cylinder,  followed  by  a  succession  of  finite-element 
atreas  analyses  with  Re  and  the  interface  pressure  being  updated  between 
each  iteration. 

An  Initial  approximation  for  Ra  for  each  cylinder  can  be  obtained  by 
performing  a  finite-element  atress  analysis  of  the  motor  with  all  in¬ 
flight  loads  present,  but  with  the  overwind  omitted.  The  derived  axial 
streaa  solution  is  then  used  together  with  Equation  6  to  obtain  the 
required  Ra  values  with  the  corresponding  interfacial  pressures  pb 
following  from  Equation  10. 

Once  the  initial  Interface  pressure  has  bean  calculated  the  isotropic 
finite-element  analysis  can  be  performed  using,  in  this  case,  all  in¬ 
flight  loads  together  with  the  calculated  Interface  pressures  applied 
externally  to  the  regions  previously  reinforced  by  the  fibre  overwind. 

An  improved  estimate  for  the  axial  stress  is  then  obtained  and  from  this 
a  modified  value  for  Ra  for  each  cylinder  is  determined  from  Equation  6. 
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The  updated  Interface  preasurea  then  follow  from  Equation  10.  The  whole 
proceaa  la  repeated  until  the  difference  between  aucceaalve  stress 
solutions  la  less  than  a  few  per  cent.  Computations  have  shown  that 
convergence  la  obtained  within  four  or  five  Iterations.  This  procedure 
la  shown  schematically  In  Figure  1. 

3.2  Fhotoelastlc  Analysis 

Nozzle  stresses  were  determined  experimentally  using  the  three-dimensional 
'frozen-stress'  photoelastic  technique  (Refs  12  and  13).  The  models 
(twice  full  size)  were  machined  from  cast  blanks  of  Araldlte  CT200  and  a 
separate  model  was  used  for  each  form  of  loading  (see  below).  The  results 
were  then  combined  by  simple  superposition  to  obtain  the  final  stress 
state.  Each  stress-frozen  model  was  sliced  In  both  the  meridional  and 
clreuatf erentlal  planes  and  fringe  readings  were  taken  In  a  standard 
diffuse-light  polarlscope  at  the  Inner  and  outer  edges  of  each  slice. 
Standard  techniques  were  used  for  material  calibration  and  the  motor 
stresses  were  obtained  by  proportional  scaling  as  outlined  below. 

The  approach  for  each  of  the  four  forme  of  loading  studied  Is  dicussed 
below. 


3.2.1  Inertial  Loading 

Stresses  due  to  the  Inertia  loadings  resulting  from  a)  the  rocket  motor's 
own  mss  and  b)  the  fins  were  studied  separately.  The  former  were  obtained 
from  a  model  subjected  to  s  simple  dead  weight  axial  loading  (see  Fig  2). 
The  Inertia  load  at  any  transverse  plane  In  the  motor  Is  equal  to  the 
product  of  the  mass  of  that  part  of  the  motor  between  the  plane  and  the 
outlet,  and  Its  acceleration.  In  order  to  obtain  the  Inertia  stress  at 
any  particular  transverse  section  of  the  motor,  the  model  stresses  were 
Increased  by  a  factor  of  four  (to  allow  for  model  scale)  and  then  by  the 
ratio  of  the  calculated  inertia  load  at  the  section  to  the  load  applied 
to  the  model. 

Stresses  due  to  the  fln-loadlng  were  obtained  from  a  model  subjected  to 
dead  weight  loads  applied  to  metal  strips  fitted  Into  each  of  the  four 
fin  slots  nearest  the  nozzle  outlet.  Motor  stresses  were  obtained  by 
proportionate  seeling. 

3.2.2  Fra s sure  Loading 

The  non-uniform  Internal  pressure  In  the  motor  was  slMilated  In  the  model 
by  a  hydrostatic  pressure  loading  using  a  combination  of  glycerine  and 
Hood's  metal  (a  low-meltlng-polnt  alloy  which  Is  liquid  at  the  stress- 
f reccing  temperature).  The  levels  of  the  two  'liquids'  Inside  and  outside 
the  model  (see  Fig  3)  were  such  that  a  two-slope  linearly  varying  pressure 
was  obtained  In  the  nozzle  region  of  the  motor  and  a  constant  pressure 
was  obtained  In  the  cylindrical  portion.  The  relative  denslcles  of  the 
liquids  were  such  that.  In  relative  terms,  this  pressure  distribution  was 
closely  similar  to  that  In  the  motor  (Pig  4).  Motor  stresses  were  obtained 
from  model  stresses  by  simple  scaling. 
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3.2.3  Overwinding 

The  pre-stress  resulting  from  the  overwinding  process  was  simulated  by 
shrink-fitting  *  series  of  Araldlte  rings  on  the  parallel  cylindrical 
section  of  the  aodel.  The  ring/cylinder  interference  was  chosen  so  as  to 
give  acceptable  fringe  orders  in  the  stress-frozen  aodel;  the  rings  were 
heated  to  135*C  to  allow  them  to  be  fitted  on  to  the  cylinder.  An 
expression  la  available  (Ref  3)  for  the  mean  circumferential  stress 
produced  In  the  cylinder  in  terms  of  the  relevant  motor  variables.  In 
obtaining  tbe  motor  stresses  from  the  model  results,  the  model  stresses 
were  scaled  by  the  ratio  of  the  calculated  mean  circumferential  stress  In 
the  motor  to  the  observed  mean  circumferential  stress  In  the  cylindrical 
portion  of  the  model. 

4  COMPARISON  OP  FINITE  ELEMENT  AND  PWOTOELASTIC  STRESS  RESULTS 


Finite-element  and  photoelastic  results  have  been  compared  for  a  typical 
motor  case  using  representative  material  property  data  and  ln-flight 
loads  (see  Tab  1).  A  typical  non-uniform  pressure  distribution  acting  on 
the  Inner  surface  of  the  motor  and  nozzle  la  shown  by  the  solid  line  in 
Figure  4. 

Because  of  the  complexities  associated  with  the  motor  geometry  and  In¬ 
flight  loadings,  certain  simplifying  assumptions  were  Introduced  Into  the 
finite-element  analysis.  The  fin  slots  and  location  holes,  which  could 
not  be  modelled  by  the  axlsymmetric  finite-element  program,  were  omitted 
and  replaced  by  solid  sections  In  the  finite-element  mesh.  Furthermore, 
to  model  the  Influence  of  the  fins  an  element  of  representative  mass  was 
positioned  adjacent  to  the  fins'  rear  location  point.  The  final  finite- 
element  mesh  employed  (see  Fig  5)  represented  a  compromise  between  solution 
accuracy  and  computer  storage.  Saint  Venant’s  principle  was  invoked  and 
the  parallel  section  of  the  motor  was  truncated  one  diameter  In  length 
from  the  end  of  the  overwound  portion  of  the  motor.  A  comparison  between 
the  finite-element  and  thin-walled  cylinder  solutions  showed  that  the 
Influence  of  the  nozzle  did  not  extend  beyond  that  point. 

The  resulta  of  the  finite-element  and  photoelastic  stress  analyses  are 
presented  in  Figures  6  to  11.  Selected  plots  of  the  circumferential  and 
meridional  stress  distributions  at  the  Inner  and  outer  surfaces  of  the 
motor  for  the  same  croas  section  are  Illustrated.  A  plot  of  the  deformed 
ehape  obtained  from  the  finite-element  analysis  Is  alto  presented  in 
Figure  12. 

In  Figures  6  and  7  the  circumferential  and  meridional  stresses  are  plotted 
respectively  for  the  overwinding  and  inertia  loading  conditions.  These 
two  quantities  represent  the  dominant  strasses  Which  era  produced  for  theae 
two  load  eaaaa.  In  Figure  6,  where  the  circumferential  stress  distributions 
obtained  from  both  enelyees  ere  plotted  for  the  Innar  eurface,  the  etress 
prediction#  ere.  In  general,  In  good  agreement.  However,  differences  ere 
observed  for  the  region  at  the  end  of  the  overwound  section  which  Is 
dose  to  the  nozzle  inlet.  In  the  finite-element  analysts  the  thin-walled 
cylinder  equations  were  need  to  obtain  the  interface  pressure  np  for  this 
region  end  consequently  'end  effect'  errors  may  have  been  Introduced. 
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This  nethod  of  evaluating  pj,  Is  an  Intrinsic  feature  of  the  approximate 
f lnlte-eleoent  procedure  and  the  introduction  of  these  errors  Is  unavoidable. 

In  figure  7,  where  Meridional  stress  distributions  at  the  Inner  surface 
of  the  rocket  notor  case  are  plotted  for  the  Inertia  loading  condition, 
it  can  be  seen  that  Che  finite-element  stress  predictions  are  again  In 
reasonable  agreement  with  the  photoelastic  results.  However,  unexpected 
differences  occurred  In  the  thln-wallad  parallel  section  of  the  motor, 
and  since  the  finite-element  and  analytical  solutions  are  In  close 
agreement  this  difference  has  probably  been  Introduced  through  a  deficiency 
in  the  Inertia  modelling  used  in  the  photoelastic  analysis. 

In  the  remaining  figures  (Figs  8-11)  both  the  meridional  and  circumferential 
stress  distributions  are  plotted  for  both  the  Inner  and  outer  surfaces 
for  the  combined  loading  case,  le  overwinding.  Inertia  and  pressurlsatlon. 
Since  the  pressurlsatlon  analysis  is  Included  In  the  combined  loading 
analysis  a  separate  assessment  of  that  loading  case  has  not  been  presented. 
Again,  there  Is  good  agreement  between  the  finite-element  and  photoelastic 
results.  However,  the  following  points  are  worth  noting: 

a.  In  Figures  8  and  9,  where  the  circumferential  stress 
distributions  are  plotted  for  the  Inner  and  outer  surfaces 
respectively,  differences  In  the  peak  stresses  are  observed. 

These  differences  are  attributable  to  the  coarseness  of  the  finite- 
element  mesh  employed  adjacent  to  sudden  changes  In  motor  geometry 
and  the  three-dimensional  effects  of  the  fin  slots  and  location 
holes  Included  in  the  photoelastic  analysis. 

b.  For  combined  loading  conditions  there  Is  also  good 
agreement  between  th*  two  stressing  techniques  for  the  prediction 
of  the  meridional  stress  distribution  at  the  Inner  (see  Fig  10) 
and  outer  (see  Fig  11)  surfaces.  The  differences  which  were 
present  in  the  Inertia  analysis  have  been  carried  through  Into 
the  combined  loading  analysis.  The  photoelastic  predictions  for 
the  thln-walled  section  consequently  differ  slightly  from  that 
given  by  the  finite-element  solution. 

c.  Although  the  Interface  pressure  between  the  first  fibre 
layer  and  the  case  can  be  evaluated  theoretically,  It  Is  nevertheless 
an  empirically  unknown  quantity.  Consequently  photoelastic  stress 
solutions  could  not  be  obtained  for  the  outer  surface  of  the 
overwound  portion  of  the  notor. 

5  CONCLUSIONS 

Stress  analysis  procedures  using  the  finite-element  and  photoelastic 
techniques  for  analysing  flbre-relnforced  metallic  rocket  notor  cases 
have  been  presented.  These  techniques  enabled  analyses  of  the  overwinding 
process.  Inertia  and  pressurlsatlon  due  to  firing  to  be  performed. 

Certain  nullifications  and  aesmsptlons  had  to  be  Introduced  Into  both 
methods.  When  the  two  procedures  were  applied  to  the  stress  analysis  of 
a  typical  overwound  rocket  motor  case  reasonable  agreement  ms  achieved 
for  all  load  cases  studied. 
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TABLE  1  Motor  Case  Material  Properties  and  In-Flight  Loads  Data 


Young's  modulus 

Fibre 

: 

93GPa 

Case 

: 

73.2GPa 

Poisson's  ratio  of  case 

: 

0.3 

Yield  stress  of  case 

: 

f'7(MPa 

Density  of  case 

: 

2700kg /m3 

Thicknesses 

Fibre 

: 

0. 1 1mm 

Case 

: 

1.95mm 

Radii 

Interface 

. 

46.9mm 

Effective 

2 

34. 6mm 

Applied  winding  strain 

: 

0.33X 

Applied  winding  strain 
ftimber  of  layers 


20 
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FIG.  1 


FIG.  1  APPROXIMATE  FINITE-ELEMENT  SOLUTION  PROCEDURE 
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30  x  21b 


<(>  200mm 
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FIG.  3 


FIG.  3  NON-UNIFORM  PRESSURE  SIMULATION  PROCEDURE 
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FIG.  4  PRESSURE  DISTRIBUTION  THROUGH  NOZZLE 
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RENTIAL  STRESS  vs  AXIAL  STATION:  OVERWOUND  STATE  (INNER  SURFACE) 
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FIG.  7 
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Fia  7  MERIDIONAL  STRESS  vs  AXIAL  STATION:  INERTIA  LOADING  (INNER  SURFACE) 
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FIG.  8 
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AXIAL  POSITION  mm 

FIG.  8  CIRCUMFERENTIAL  STRESS  vs  AXIAL  STATION:  COMBINED  LOADING  (INNER  SURFACE) 
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FIG.  9 
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AXIAL  POSITION  mm 

FfG.  9  CIRCUMFERENTIAL  STRESS  vs  AXIAL  STATION:  COMBINED  LOADING  (OUTER  SURFACE) 
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FIG.  12 
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ANNEX  A 


INTERFACE  PRESSURE  CALCULATION 


i 


The  Interface  pressure  (pjj)  between  the  first  fibre  layer  and  the  case 
consists  of  two  coiqponents,  one  due  to  winding  (p^)  and  one  due  to  the 
application  of  the  Internal  pressure  (p^p)* 

The  winding-induced  component  p^u  can  be  readily  obtained  from  the  pre¬ 
stress  produced  during  winding  through  the  thin-walled  cylinder 

equations  in  the  form 


*w 


(Al) 


which,  when  using  Equation  2  to  eliminate 


°ec 


(n) 

» 


modifies  to: 


Ec  i  ^ n 
Pbw  * - R - 


(A2> 


In  principle  it  would  also  seem  admissible  to  use  the  thin-walled  cylinder 
approximations  to  calculate  the  pressure-induced  Interface  pressure 
conponent  pj,p  to  give 


Pbp 


09cp(n)tc 
P  “  Pbw  "  — i - 


(A3) 


and  hence 


Pb  "  Pbw  +  Pbp 

(A4) 

However,  the  thin-walled  cylinder  equations  are  approximations  of  the 
exact  solution,  as  given  by  Lame  (Ref  14).  Although  the  error  in  pj,p 
is  consistent  with  these  approximations,  typically  less  than  a  few  per 
cent,  it  has  been  found  that  under  certain  conditions  considerable  errors 
In  the  circumferential  stress  calculations  can  occur  when  this  estimated 
value  for  pt,p  is  used.  If  an  error  analysis  is  performed  for  the  circum¬ 
ferential  stress  og  to  assess  the  effects  of  small  errors  in 
p~,,  on  a g  (see  Annex  B)  the  following  relationship  can  be  derived 

A09  _  APbp 
OQ  Pbp”P+Et  c 

R  (A5) 


«  where  A09  and  Apt,p  are  the  respective  errors  in  09  and  p^p.  For  a  thin- 

walled  cylinder  clrcisaf  erentially  reinforced  with  a  fibre  overwind  the 


Al 
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interface  pressure  pj,p  can  approach  the  internal  pressure  p,  resulting 
in  the  equation 


Age  .  **Pbp 

a8  tpPbp 

(A6) 

From  the  above  equation  it  is  evident  that  vhen  p  and  pbp  are  approximately 
equal  the  relative  error  in  og  is  proportional  to  the  relative  error  in 
Pbp  with  the  constant  of  proportionality  given  by  R/tc .  Since  for  thxn- 
walled  cylinders  R/tc  is  typically  greater  than  ten  it  follows  that  small 
errors  in  the  interface  pressure  predictions  can  induce  large  errors  in 
the  circumferential  stress  calculations. 

If  the  Interface  pressure  is  significantly  different  from  the  Internal 
pressure,  as  is  the  situation  during  winding,  the  relative  errors  in  the 
interface  pressure  predictions  are  comparable  to  those  of  the  circumferential 
stress  calculations.  It  is  for  this  reason  that  it  is  admissible  to  use 
the  thin-walled  cylinder  approximations  to  calculate  p^w.  la  Equation  A2. 


In  view  of  the  previous  discussion  it  is  essential  therefore  that  the 
pressure- induced  interface  pressure  pbp  is  evaluated  using  the  Lame 
equations.  Omitting  details  it  can  be  easily  shown  from  a  compound 
cylinder  analysis  that  the  required  Interface  pressure  pj,p  is  obtained 
from 


2V 


Ec<V-Ri2  ) 


Pbp 


r  Rh2+R<2  +  R„2+Rh2  1 


(A7) 


where  Rj  and  are  the  internal  and  interface  radii  respectively  of  the 
metallic  case  and  RQ  is  the  external  radius  of  the  fibre-reinforced  metallic 
rocket  motor  (ie  R„»Rb+ntf). 
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hand  If  pfep  equals  p  Chen  Che  reladve  error  In  a 3  becones  proportional 
Co  the  relative  error  in  pj,-  with  the  constant  of  proportionality  given 
by  R/tc.  Since  for  thin-walled  cylinders  R/tc  will  always  be  significantly 
greater  than  unity  the  relative  error  in  03  for  a  small  error  in  pbp  can 
become  large. 
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